SO(5)-Symmetric Description of the Low Energy Sector of a Ladder System 
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We study a system of two Tomonaga-Luttinger models coupled by a small transverse hopping (a 
two-chain ladder). We use Abelian and non-Abelian bosonisation to show that the strong coupling 
regime at low energies can be described by an SO(5)i WZW model (or equivalently 5 massless 
on : Majorana fermions) deformed by symmetry breaking terms that nonetheless leave the theory critical 

at T = 0. The SO (5) currents of the theory comprise the charge and spin currents and linear 
combinations of the so-called pi operators (S.C. Zhang, Science 275, 1089 (1997)) which are local in 
terms both of the original fermions and those of the effective theory. Using bosonisation we obtain 
the asymptotic behaviour of all correlation functions. We find that the 5 component "superspin" 
£f>^ , vector has power law correlations at T — 0; other fermion bilinears have exponentially decaying 

correlations and the corresponding tendencies are suppressed. Conformal field theory also allows us 
to obtain the energies, quantum numbers, and degeneracies of the low-lying states and fit them into 
. deformed SO (5) multiplets. 

C/3 



I. INTRODUCTION 



One of the most characteristic features of the high-T c Cuprates is the proximity of antiferromagnetic (AF) and 
superconducting (SC) phases as a function of doping. As a result, much of the theoretical effort has focused on trying 
to consistently treat the insulating-underdoped-optimally doped region of the phase diagram, in which AF and SC 
tendencies compete and may have strong fluctuations. 

An interesting recent proposal is that of Zhang .El He suggests that the simplest way of unifying AF and SC in the 
Cuprates is to introduce a new five-component vector order parameter consisting of the three component staggered 
magnetisation, and two components associated with the real and imaginary parts of the d-wave SC order parameter. 
Clearly this new concept is only useful if there exists some kind of symmetry (higher than the known SO(3)®U(l)) 
which relates the AF and SC sectors. His suggestion is that an approximate SO(5) symmetry emerges in the low 
I . energy sector (SO (5) because the new composite order parameter has five components and transforms like a vector). 
[ — ' If true, this would allow the construction of an SO(5) quantum nonlinear er-model to explain the low-energy dynamics 
ON . of the high T c materials. This could explain the form of the phase [diagram, and the so-called 7r-modcB 

However, there have been several criticisms of this theory. Sometl have focused more on the details of microscopic 
calculations in the framework of the t-J or Hubbard models. Others have added several physical objections.Q One 
response to these criticisms has been to attempt to construct concrete examples of extended microscopic Hamiltonians 
which manifestly have an SO (5) symmetryia But knowing the Hamiltonian does not necessarily tell us much about 
the low energy behaviour. 

In this paper we study a two-chain ladder Hamiltonian that is related to popular two-dimensional models of the 
Cuprates. One of the reasons that ladder systems have attractecLsuch attention is that many experimental realisations, 
of these systems are very closely related to the high T c materials^, and some have even exhibited superconductivityfl'El 
From a theoretical point of view, powerful non-perturbative techniques such as bosonisation and conformal field theory 
(CFT) exist in one dimension. This offers hope of starting with a microscopic Hamiltonian and ending up with a 
tractable effective field theory. In this paper it is not our purpose to comment on the general validity of the SO (5) 
. . . i idea but to explicitly study a simplified and more tractable model. n .— . 

There is a large body of literature on two-chain and ladder systemstrEZl (for a review see Ref. ||) . Using a combination 
of weak coupling RG and bosonisation, the phase diagram has been intensively investigated. These analyses reveal 
that for small interchain hopping there are interesting strong coupling phases. However, whilst Abelian bosonisation 
and weak coupling RG are good for determining the phase diagram, they do not explicitly respect the symmetries of 
the system, nor do they provide detailed information about the correlations. In this paper we explore in more detail 
the strong coupling region of a two chain ladder system, taking care to preserve the full non-Abelian symmetries and 
obtain the correlations. 

It is well known that many 2 chain ladder systems are spin liquids; that is, they exhibit a spin gap for a wide range 
of different fillings and couplings. This is because the Luttinger liquid is a quantum critical system, and as such, 
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highly unstable to perturbations such as interchain coupling. In general there are a number of relevant couplings 
which can drive the system into a spin gap phase (an explicit example is discussed in section VI). However, in this 
paper we study a simplified system in which there is no backscattering and as a result, no spin gap. This model is of 
interest because it displays remarkable similarities to some aspects of the Zhang proposal in 2 dimensions .EJ 

The model we consider is a system of 2 spinful Tomonaga-Luttingcr (TL) models in the repulsive regime, coupled 
by a small interchain hopping. This corresponds to the case of no backscattering and was studied in Refs. [If], [Tt]. 



We demonstrate that the hopping only generates couplings in a certain sector of the theory (which we call "flavour"), 
freezing it out of the effective action at energy scales below tj_. In agreement with the above references, we find that 
this leaves a critical (at T=0) spin and charge sector with conformal charge 5/2. However, we go on to show that 
this can be represented as a system of 5 massless Majorana fermions, or equivalently, an SO(5)i Wess-Zumino-Witten 
(WZW) model, deformed away from the symmetric point by marginal current-current interactions. These SO (5) 
breaking terms are associated with spin charge separation (spin and charge velocities not equal v s ^ v c ) and the 
anomalous charge exponent (K c ^ 1), which distinguish the spin and charge sectors. Thus the system is never exactly 
SO (5) symmetric except in the trivial noninteracting case. Nonetheless, this representation does have strong analogies 
with the Zhang proposal in 2d; the physics can be understood using an SO (5) symmetric er-model with symmetry 
breaking terms. In this way we obtain the asymptotic behaviour of all correlation functions; the correlations of the 5 
component "superspin" are enhanced (power law at T=0); we obtain their scaling dimensions. Other fermion bilinears 
die away exponentially fast. 

Sections II- V are concerned with an analysis of this model, including its detailed symmetric description, the relevant 
currents, the 7r-operators, its correlations and low lying multiplets in the excitation spectrum. One important way in 
which the system we are studying differs from that considered by Zhang is that we are away from half-filling, which is 
a very special point in Id— Exactly at half-filling it is necessary to consider the Umklapp term, which causes a Mott 
gap in the charge sector.lla Then the low energy effective Hamiltonian is simply a pure spin Heisenberg model (with 
exchange J ~ 4t 2 /U in the case of the repulsive Hubbard model at strong U). We comment further on this difference 
in section IV. 

In section VI we finally consider the case of two coupled Luttinger liquids, which differs from the previous model in 
that it includes marginal backscattering terms. An example of this is provided by some regions of the phase diagram 
of a system of 2 Hubbard chains coupled by single particle hopping. In this more physical case, we show in detail how 
the additional marginal terms cause a spin gap to appear in agreement with Refs. ^ [H], and numerical work such 
as Ref. p(| Then the spectrum and correlations are as in Ref. ^l]; there is a spin gap but the charge sector remains 
gapless. 

Finally, we conclude. There is also an appendix which sketches out a bosonisation prescription that enables us to 
calculate the correlation functions of fermion bilinears. 



Many systems of interacting one-dimensional fermions away from half filling fall into the Luttinger liquid universality 
class. That is, they exhibit spin-charge separation, gapless excitations, anomalous power law correlations and the 
absence of a quasiparticle pole (see Ref. |22| for a recent review, and references therein). For example, the one- 
dimcnsional repulsive Hubbard model away from half-filling is known from its exact solution to be a Luttinger liquid 
all the way from U = to U = oo, as is the t-J model for small enough J/t. 

One of the simplest two-chain models of this type that can be written down consists of two Tomonaga-Luttingcr 
(TL) models (labelled by a chain index i = 1,2) coupled by a small transverse hopping t± -C t: 




II. A SIMPLE MODEL 



H = tf TL (l) + H Th {2) + H± 



(1) 



where the TL Hamiltonian is a sum of three pieces (Hq + if 2 + H4): 






(2) 



The current (or density) is simply defined as 
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3a,i — Ra, t Ra,i 3a,i ~ ^a,i^">i (3) 

and the electrons fields R a ,i and L a ^ are slowly varying on an atomic scale: the electron annihilation operator at site 
x, chain i and spin a may be expressed as 

c a ,i{x) = R a ,i{x)e ikFX + L a ,i{x)e- ikFX (4) 

In terms of these fields, the simple interchain hopping term becomes 

H±=t± f dx (rI,i{x)R<*,2{ x ) + L i,x{x)L a<2 (x) + h.c.) (5) 

For simplicity, we have assumed that the Hamiltonian is invariant under spin rotation, and so the coupling constants 
52 and (74 are the same for parallel and antiparallel spin configurations. Normal ordering is assumed throughout in 
products of local fields (definition of currents, Hamiltonians, etc.). 

It is worth making a quick observation about the difference between the terms TL liquid and TL model: The TL 
model is an idealised and specific Hamiltonian, written down in Eq. (||). It has a perfectly linear dispersion, an infinitely 
deep Fermi sea, has only density-density interactions i^-pd is exactly solvable for all values of the coupling constants 
(the model is unstable beyond a critical value of c^)-^ The TL liquid (which is the generic state corresponding to 
many realistic Hamiltonians like the Hubbard model away from half filling) differs in that the dispersion is no longer 
exactly linear, and the Fermi sea no longer infinitely deep. But from our point of view the most important difference 
in the low energy sector is the presence of marginally irrelevant couplings (backscattering) . In a single chain system 
these are not very important when repulsive - they simply give logarithmic corrections to the correlation functions. 
In section VI we will study the effect of these additional terms in the two chain system, in order to establish the 
behaviour of the more realistic coupled TL liquids, but for the moment, we will restrict our attention to the simpler 
case of coupled TL models. 

The model (|), even though it is made of TL models, is not exactly solvable because of the interchain hopping. 
However, we will argue presently that the model segregates into three different sectors, respectively associated with 
charge, spin, and "flavour" , and that the combined effect of interchain hopping and interactions is to make the flavour 
sector massive, leaving only the charge and spin sectors critical (i.e., gapless). To each sector one may associate 
current operators, expressed as bilinears of the electron fields: 

charge: J R (x)=Y j W a ^(x)R ati {x) 

spin: 3 R (x)= - ^ R l,i( x ^ (Ta R l 3 -^ 

i,at,/3 

flavour: Ir(x)= - ^ R^ i {x)cn j R a ,j{x) (6) 



where er is the vector of Pauli matrices (left-moving currents are defined similarly) . These currents have the following 
commutation relations (they may be derived from Wick's theorem): 

[J R (x),J R (y)] = --S'(x-y) 

7T 

[Jr( x \ J h R (y)} = -^ ab S'(x -y) + is abc J R (y)5(x - y) 

[IrW, 4(2/)] = ~^ 6ij6 '( x ~V) + m jk I R {y)5{x - y) (7) 

and currents of different types (i.e., charge, spin, and flavour) commute. Thus in the language of non-Abelian boson- 
isation, tter-eharge current obeys a U(l) Kac-Moody algebra, and the spin and flavour currents obey SU(2)2 =SO(3)i 
algebras. ErEZl It is simple to show that the Hamiltonian (Q) may be expressed as H = Hq + V c + Vf , where cynly the 
above currents appear. This is just a matter of taking careful account of point-splitting and normal orderingO 

H = ^Jdx(j R + 3l + ll+ [R^L]) 
V c = ^Jdx {g 2 J R JL + g4(J 2 R + J 2 L )} 

V f = 2jdx {g 2 F R Il + 9i [(I R f + (II) 2 } + t ± (I R + ID) (8) 
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Therefore the model ([!]) decouples into three independent sectors (charge, spin, flavour). The important point is that 
the hopping term only involves the flavour sector, which is decoupled from the other two. The effect of interactions 
(172 and #4) on the charge sector will be a velocity renormalisation and anomalous scaling exponents (K c 7^ 1). The 
combined effect of interactions and transverse hopping on the flavour sector is more dramatic. The RG analysis of 
Rcf. ^ shows unambiguously that in the repulsive regime (K c < 1 ), the system scales to strong coupling at energies 

< t± (in the notation of Rcf. [I?] our model corresponds to initial conditions of — 0, g^ = — = go for 
i = 0, 7r, /, t, b). The combination of the small hopping term t±_ and the interaction terms leads to the generation of 
important couplings in the RG process, giving a gap in some channels. What our analysis tells us is that all of this 
physics is only happening in the flavour sector, and thus it is this sector that becomes gapped, while the total spin 
and total charge sectors remain untouched and critical. So at low enough energies the flavour sector is frozen out of 
the effective theory, and our task is simply to understand the remaining charge and spin degrees of freedom. 



III. SPINOR AND VECTOR DESCRIPTIONS 



Each electron field R a ^ or L a i carries charge, spin and flavour. The separation of the model into charge, spin and 
flavour sectors is therefore difficult to describe in terms of these operators. However, one may introduce a different set 
of Fermi fields in terms of which this separation is much more natural. To this end, we must use some representation 
theory of Lie groups. 

Let us first consider the model (|f|), but without interactions or interchain hopping (i.e., two free, decoupled chains). 
This model has SO(8) symmetry, and this may be shown has follows. Each complex field R,L may be written in 
terms of its real and imaginary parts: i? Q ^ = i?i jQ ,,j + iRi,a,i an d then, except for a total derivative, the Hamiltonian 
Hq takes the form 

H = iv F ^2 dx {R^d x R^ - L M 9 X L M ) (9) 

where the composite index fi, running from 1 to 8, stands for spin, chain and real/imaginary part. The eight Fermi 
fields i? M (or L^) can undergo an internal SO(8) rotation that leaves Ho invariant. Hence the model has a chiral 
SO (8) symmetry, ft is well known that a collection ofJV real free fermions like this is equivalent to a special kind of 
conformal field theory: a level-1 SO(N) WZW modelo Chiral SO(8) currents may be defined in terms of those real 
fermions as follows: 

1 8 

ii.y— 1 

where S™, is a matrix representation of the generators of SO(8) (A runs from f to ^N(N — 1 ) = 28, the number of 
generators). Left-moving currents are defined similarly. The charge, spin and flavour currents (||) are special cases of 
the above and correspond to specific values of the index A if the generators S£- are chosen judiciously. 

The currents (|To| ) are bilinears in the electrons fields i? M (a = 1, ... ,8). However, the SO(8)i WZW model contains 
other fields, belonging to a different representation of SO (8), in terms of which these currents are also bilinears. Among 
all SO(N) groups, SO(8) is peculiar in that its vector representation, of dimension 8, has properties identical to its 
spinor and conjugate spinor representations (also of dimension 8). Indeed, which one is called 'vector' is a matter of 
convention, dictated by the way the SO (8) symmetry breaks down to smaller SO(N) components. In order to decide to 
which SO (8) representation the electron fields belong, one must study in detail how each representation breaks down 
when the symmetry is reduced. Let us consider a two-stage symmetry breaking, in which the flavour sector, with its 
SU(2), is first segregated, and then the charge U(l) and spin SU(2) (note that U(l)~SO(2) and SU(2)~SO(3)): 

50(8) -> 50(5) ® 50(3) fl - -► 50(2) c ® 50(3) sp - <g> 50(3) fl ' (II) 

We stress that the goal of the present analysis is to fit the fields and states of the model into symmetry multiplets, 
without demanding the symmetry to be exact. In the first stage of this breakdown, the vector and spinor represen- 
tations of SO (8) are decomposed as follows (irreducible representations will be commonly denoted by bold numbers 
giving their dimensions, with an occasional superscript distinguishing between vectors (v) and spinors (s)): 

8 V -» (5,1) (1,3) 

8 s (4, 2) (12) 
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Q 



Q 




FIG. 1. The S z — Q diagrams associated with the lowest nontrivial SO(5) multiplets. 



(here, for instance, the notation (4,2) stands for a tensor product of the 4-dimensional representation of SO(5) 
with a doublet of SU(2) fl ). Since SO(5) representations are not all that familiar, we provide a pictorial view of 
the lowest nontrivial ones on Fig. [j]. The multiplet 4 is the spinor representation of SO (5), while 5 is the vector 
representation and 10 the adjoint representation, i.e., the representation of the SO(5) symmetry currents or generators. 
The decomposition of these SO (5) representations in terms of spin multiplets and charge quantum numbers is best 
appreciated on Fig. |l|. For instance, the SO(5) spinor 4 breaks down into two spin-i doublets, one with charge +1 
and the other with charge -1. On the other hand, the vector representation breaks down into a spin-1 triplet of charge 
zero and two singlets of charges ±2. 

We may now ascertain that the electron fields belong to the spinor representation of SO (8). Indeed, the lowest 
excited states of Hq, obtained by acting on the vacuum with the lowest electron creation operators, form a multiplet 
of 4 states of charge +1 and four states of charge —1. This is precisely the charge content of the spinor multiplet 
8 s , since the spinor 4 of SO (5) contains two states of charge +1 and two of charge —1, and appears twice in the 



decomposition (12), because of the flavour doublet. 

A different set of real fermions, denoted & (i — 1, ... ,8), belongs to the vector representation of SO(8). These new 
fermions are related in a complicated, nonlocal way to the original fermions. The transformation relating them may 
be explicitly obtained via Abelian bosonisation, if one takes care to preserve the anticommutation factors, but this is 
not a particularly illuminating procedure. The important point is that they are just a different basis or representation 
for the same system. These fermions obey the usual anticommutation relations , £j (y)} = SijS(x — y). The 

SO (8) currents ([h]) may also be expressed as bilinears of these fermions, albeit with the help of a different set of 
SO (8) matrices: 



J. 



A 



A characteristic feature of the vector representation is its particularly simple decomposition into charge, spin and 
flavour components: in the first stage of the breakdown (p"T|), the vector representation decomposes as 8 V — > (5, 1) © 
(1,3). In the second stage, the SO(5) vector decomposes as 5 — > (3,1) (1,2) (this time, doublets on the r.h.s. 
correspond to spin and charge multiplets, respectively). We may thus distinguish three Majorana fermions (£f, 
i = 1, 2, 3) for spin, three others (£/, i = 1,2, 3) for flavour and the remaining two (£?, i = 1,2) for charge. The spin, 
flavour and charge current then have the following expressions: 

i 

In (14) 

Jr = -iJ k $% = (15) 

It is interesting to note that these currents (and eventual SO(5) currents) are local in terms both of the electron fields 
and in terms of the above Majorana fermions, even though the fcrmion operators themselves are nonlocally related. 
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The Majorana fermions are nonetheless legitimate operators of the theory. For instance, at half-filling, if all sectors 
became gapped, the three spin fermions £| would describe the triplet of spin excitations characteristic of a gapped 
spin-1 chainBa 



IV. THE SO(5) CURRENTS AND ZHANG'S IT OPERATORS 

We have seen above that the spin and charge degrees of freedom, which make up the critical sector of the theory 
(|l|), may be described by the five Majorana fermions ££ 2 an d £1,2,3- Except for the interaction V c of Eq. (|^), the 
low-energy sector is equivalent to a level-I SO(5) WZW model with conformal charge c = |.fcj Indeed, the above 
Majorana fermions may be arranged in the following suggestive sequence: 



a = & 



6 = fi 6 = £2 U =€3 £5 = £1 



(16) 



plus corresponding left-moving fields. That the non-interacting part of the spin-charge sector is equivalent to a level-T 
SO (5) WZW model means that this part of the Hamiltonian may be simply expressed asE-3 



Ho = iv F ^2 

3 = 1 



(17) 



(here £j- denote the left-moving fields). 

It is then useful and instructive to introduce the full ten SO(5) currents. Four of those currents are provided by 
the charge and spin currents of Eq. (^). The remaining six, corresponding to Zhang's 7r-operators,EJ may be expressed 
in the continuum limit either in terms of the Majorana fermions 2 an d £1 2 3j or directely in terms of the electron 
fields. PI 

It is interesting at this point to go back to the lattice definition of the tt operators :tl 



t - V 



IE = 



,(k)ct(k+Q) (r.osUcU-k) 



E ffm,ne iQ - m ct a (m)(a a a 2 ) a/3 4(n) 

m,n,a,/3 



(18) 



where m and n are vectorial site indices (in-chain and chain index). On a square lattice at half-filling Zhang takes 
Q = (tt, 71"). On a two-chain system, away from half filling, there are two possibilities: Q = (2fc^,7r) for right movers 
and Q = (— 2kp,Tr) for left movers. 

The structure factor g(k) = cosk x — cosk y has the local form: 



+2 if (m,n) are NN on the same chain 
—2 if (m,n) are NN on opposite chains 
otherwise 



(19) 



Defining a "staggered 7r-density" : 



lit = 8 / dx e 2lkFX irl 



(20) 



we find, with the help of Eq. (0), the following expressions: 



7T+ = 



rt I,l- n 'i,2 
Jl J.,l- n 'i,2 



(21) 



Interestingly, this continuum expression for the tt operators is quite robust and does not depend too closely on the 
microscopic definition (|l8|). One might have alternatively chosen the Henley-Kohno form g(k) =sgn(cosfc;c — cosk y ) 
and this would not have changed the results, apart from derivative terms which are irrelevant in the RG sense - 
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essentially because in a two-chain ladder system there are only two k y values, or tt. Another seemingly different 
microscopic expression for the tt operators is used in Ref. |l^, but again, we have verified that the same expression 
( pl| ) is obtained in the continuum limit. 
We define the matrix l ab (x), analogously to Zhang: 

/ \ 

(22) 



4 +7T X 







lA + TTy 

7T] + TT Z 


J R 





J R 


— J x 

J R 


Jr 


- TT X ) 


-i(4 - 





V J R -i{-ni-n x ) -i(nl - 7T y ) -i(?r| - tt 2 ) Oj 

(the matrix is antisymmetric and so we only wrote down the lower diagonal). Using Wick's theorem for the electron 
fields R ay i and L a ^, we find that the l ab obey an SO(5)i Kac-Moody algebra, different from the standard SO(5) 
algebra by a quantum anomaly coming from the necessity for normal ordering with respect to the vacuum: 

[l a \x), l cd {y)] = 5(x - y)(S ac l bd (x) - 5 ad l bc (x) - 5 bc l ad {x) + 5 bd l ac (x)) + ^-S'(x - y)(S ac S bd - 8 ad S bc ) (23) 

A similar procedure gives the corresponding relationship for left-moving currents. 

How can the SO (5) symmetry currents be expressed in the Majorana (vector) representation £? A vector represen- 
tation of the SO(5) generators may be easily written down: 

t% b) =i(S°S b -S*5 b ) (24) 

The currents ( p2| ) can then be represented as follows: 

n*) = \ E ^ b) ^ ( 25 ) 

z,j=1..5 

where the 5 fermions & are numbered as in Eq. ( |l6| ) : Note that the tt operators correspond to bilinears involving one 
fermion £ from the spin sector and one from the charge sector, which again fits with the physics since we know that 
they create objects with both spin and charge. 

In the low-energy sector of the model (|l|), the spin-charge sector can be represented by an SO(5)i WZW model, 
perturbed away from the perfectly symmetric point by current-current interactions (the only interactions present in 
the spin and charge sectors in our model). So Zhang's idea of using an SO (5) er-model representation with symmetry 
breaking interactionscl is explicitly seen to be valid for this model, and the Hamiltonian for the spin and charge sectors 
can be written in terms of the SO (5) currents in the Sugawara form, analogous to the form proposed in Ref. [I]: 

H cs = H 0s + Hoc + V c 

= ? /^E{( rb ) 2 + ( ra ") 2 }+ /^{^ 15 F 5 + 54 ((/ 15 ) 2 + (P) 2 )} (26) 

^ a<b J 

One notable difference between this system and that considered by Zhang is that here we are working away from 
half-filling; the n operators are defined slightly differently, to carry momentum (±2fcF,7r). The chemical potential 
term in the Hamiltonian, which in 2d breaks SO (5), here (due to perfect nesting) simply renormalises the wavevector 
fcj?; momentum is still conserved and the algebra ( p3| ) still closes because operators carrying 2k p only give non-zero 
expectation values when combined with operators carrying — 2kp. If we were to work at half-filling, there would be 
an additional Umklapp scattering which would lead to a Mott charge gap.Ea 

Sinm,there are tt operators in this system, one may also ask whether there is a well-defined 7r-resonance as claimed 
in 2dBEJ If this is so, the commutator of the Hamiltonian with the n operator will be proportional to the n operator .El 
It can easily be seen that this will not be the case. Later on, in section V, we obtain bosonised forms for these 
operators. Then one can see that in Fourier space, the correlator of tt operators does not have a simple pole; their 
effect is not to generate a single well-defined triplet excitation but a shower of unconfined spinons and holons. 



V. BOSONISATION 



The low-energy sector of the model (|1|), i.e., the perturbed SO (5) WZW model of Eq. (26) is exactly solvable, in 
the sense that we may find the exact energy levels of low-lying states and the long-distance correlations of various 
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operators. We will first indicate the physical content of the SO (5) WZW model without perturbations, and then see 
explicitly how the interactions 172 and 34 separate spin and charge sectors, affect correlation exponents and deform 
low-energy SO(5) multiplets. „ 

In the language of conformal field theory,t3 particularly useful when dealing with critical theories, each WZW model 
contains a finite number of primary fields, having well-defined conformal dimensions A and A. An operator A with 
such conformal dimensions has the following dynamical correlations: 



1 

- vt) 2K {x + vty 



(A(x, t)A(Q, 0)) ~ - _^ 9/w _ ; (27) 



The level-1 SO(5) WZW model has two primary fields: a five-component vector field £ of conformal dimension 
A = i and a four-component spinor field h of conformal dimension jg. Under SO (5) rotations, these fields transform 
respectively in the vector and spinor representations of SO(5). Of course, the field £ is made of the five Majorana 
fermions (|l6|), whereas the field h is what is left of the original electron fields i?^, originally in a spinor representation 
of SO(8), after the flavour sector has been gapped out. Freezing out the flavour part has had the effect of decreasing 
the conformal dimension of the spinor field from \ (in SO (8)) to ^ (in SO (5)), thus making it more relevant. In 
addition to these primary fields, the ten SO (5) currents (|22| ) also play a crucial role in the theory and their correlations 
may also be exactly calculated. 



A. Spin-charge separation 

When the interactions of ( p6| ) are turned on, the SO (5) symmetry is explicitly broken and the spin and charge 
sectors of the theory separate. The spin sector, unaffected by the interactions, becomes a level-1 SO(3) WZW model, 
which is the same as a level-2 SU(2) WZW model. The SU(2)2 WZW theory contains two primary fields: a spin 
triplet (or vector) £? (i = 1, 2, 3) with conformal dimension |, and a spin-i (or spinor) field g a (a =|, J,), of conformal 
dimension ™ . Products of the left- and right-moving parts of g are commonly arranged in a 2 x 2 matrix: 



q = I 919J 5T51 
\9l9l 9191 



(28) 



The charge sector becomes a U(l) theory, which may be described by a single boson field <I> C . The effect of the 
interactions on the charge boson is simply to change the spectrum, p£ anomalous dimensions (K c ^ 1) and the theory 
remains critical. Since Abelian bosonisation is fairly standard,E3l£3 we shall only state a few results. The charge 
boson $ c may be written as the sum of right and left parts: $ c = <p c + <f> c . Defining the dual field 9 C as 9 C — <j) c — cj) c , 
the charge Hamiltonian may be written as 



H c = — 
2 



dx 



K c (d x 9 c ) 2 + — {d x $ c ) 2 



(29) 



where 



K r = 



irv F + g 4 - g 2 
irv F + g 4 + g 2 



^+f) 2 -(?) 



(30) 



For more general lattice Hamiltonians that are Luttinger liquids in the low energy sector, the parameters v c and K c 
depend in a more complicated way upon the original couplings, so in the following analysis, one can just treat them 
as independent parameters whose precise value depends upon the original model. 

If <72 = 0, there are no anomalous exponents (K c — 1) and the scaling fields e ±lv/ ^' <: represent right-moving 
fermions of conformal dimensions (^,0). When 172 is turned on, the original charge current (I 1 5 = Jr) is no longer 
conserved but becomes a linear combination of currents that are still conserved (see eg. Rcf. p2| ): 



Jr = 3R cosh § - j L sinh d (K c = e 



(31) 



where jr and jh have respectively conformal dimensions (1,0) and (0,1). The chiral components <j) c and <p c mix 
through the same Bogoliubov transformation and the original fermion operators e ±J v / 47r0c aC q U j re a left conformal 
dimension: 
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A = I (Kc + 1/K C ) + \ A = i (K c + 1/K e ) - \ (32) 

Although the above results are very well known, it is worth pausing over them for a moment. They show that the 
interaction strengths g2 and <?4 are not relevant energy scales in the low energy theory. They only appear as ratios 
g2,i/vF in the renormalisation of the velocity v c and the anomalous exponent K c . This is a very non-perturbative 
result. If we recall the exact solution of the one-dimensional Hubbard model away from half-filling, we know that 
it is a Luttinger liquid for all U > from to ooH3 In this range, K c varies from its noninteracting value of 1, to 
1/2 at U = 00. Even when the on-site repulsion is infinite, its effect in the low energy sector is just a fairly small 
renormalisation of the anomalous exponent! The theory is still critical with gapless spin and charge excitations. 

Some critics of the Zhang SO (5) proposal have claimed thatbecause of the strong on-site repulsion in the Hubbard 
model, the 7r operators cannot create low energy excitationsu The argument is essentially that one is forced to put 
two electrons on the same site, which costs an energy of order U . The reason that the criticismS may be too simplistic 
is first of all that it is a single-particle argument, whereas the low energy excitations of this system are many-body 
collective phenomena, and secondly that it is a short length scale argument which may have some validity in the 
U.V.; but we are interested in the low energy I.R. behaviour which is quite different in a non- Fermi liquid such as 
the TL liquid. Even if much of the spectral weight is shifted to high energies there is still some at low energies and 
this is what dominates the low energy theory. Given that the two-dimensional Cuprates arc examples of non-Fermi 
liquids, it cannot be ruled out a priori using these arguments that even in the presence of strong on-site repulsion, 
the 7r-operators may generate low energy excitations (at least when one is slightly away from half- filling). 

Let us then see what happens to the SO(5) currents and primary fields after spin-charge separation. Three of the 
five components of £ become a spin triplet (£ s ) and the remaining two are simply cos(-\/47rV> c ) and sin(\/47r<^c)- Out 
of the 10 SO (5) currents, six - the 7r operators - are no longer conserved currents and may then be expressed as 
products of SU(2)2 P ' fields with charge fields. Schematically, 



Tr.Ti-t ~e ±iVSf ^ (z) (g)£ s ^) (33) 



When K c 7^ 1, the conformal dimensions of the 7r operators are no longer (1,0), but rather, from Eq. (|32j) and since 
the field £ s (z) has conformal dimensions (|,0), 

A = I {Kc + 1/K C ) + 1 A = i (K c + l/K c ) \ (34) 

Thus, the 7r operators are no longer conserved currents, as expected. 

As mentioned above, the spinor representation 4 of SO(5) factorizes into a pair of SU(2) doublets of charges ±1 
when SO (5) is broken. The spinor field h may thus be factorized as 

m (™tf£0 (35) 

where g is the SU(2) spinor mentioned above and the boson factors have conformal dimension \. The decomposition 
described here can be rigorously proven by checking the corresponding commutators with the currents. We obtained 
it differently, by the method of affine characters (see Ref. |2^), which we will not explain in detail here, since the 
coincidence of conformal dimensions and components is sufficiently convincing for our purpose. 



B. The SO(5) order parameter 

One of the interesting operators of the SO(5) WZW modeL(and of its perturbed version) is a continuum version 
of Zhang's five-component order parameter n a (a = 1, . . . , 5) .EJ This operator can be defined in terms of the original 
electron fields. The components ".2,3,4 correspond to the staggered magnetisation and the components nx t s to the 
d-wave superconducting order parameter. The staggered magnetisation is defined as 

nq= 4(le + Q)ff«/jc/j(k) (36) 

k,a,/3 

Picking Q = (2fci?,7r) and using Eq. (^), we find 

k,n,P 







The Q = (—2kp, ir) component of the magnetisation is just the Hermitian conjugate of the above. The d-wave order 
parameter D is defined as 

/ dxtf =Y^ 9m, n c\ (m)cj (n) (38) 
where g m ,n has been introduced in Eq. (|l^). Taking the continuum limit, we find 

vt = E ( R \A,i + L \,A,i) - R \,A,2 - L \,A,2 - r \Aa - L \A ( 39 ) 

i=l,2 

The combinations T> + and i(T> — T>^) then correspond to n\ and 715, respectively. 

An expression of the order parameter n a in terms of the scaling fields h or £ would be more useful, since the flavour 
part is not explicitly absent from the above. Such an expression is difficult to obtain in a systematic way from the 
above expressions; but one can infer what it has to be (this result can be confirmed by Abelian bosonisation; see the 
appendix). Clearly, n a should be a bilinear in h or £, with equal left and right conformal dimensions. Let us consider 
the following SO (5) tensor products: 

4 4= 105010 

5(g)5 = l©10©14 (40) 

This means that a bilinear in £ (5 components) cannot transform as a vector of SO (5), whereas a bilinear in h (4 
components) can. We thus seek an order parameter of the form 

n a = T^h~h J (41) 

where the five 4x4 matrices T a must transform as a vector of SO (5) when h and h are acted upon by a 4 x 4 unitary 
representation of SO (5). If we denote by i ah a 4 x 4 representation of the SO (5) generators, this requirement amounts 
to 

[t b , r c ] = i(s ac r b ~ s bc r a ) (42) 

Experience with the Lorentz group and Dirac matrices may guide us here. If a set of five matrices T a obey the Clifford 
algebra {r a ,r 6 } = 2S ab , then it is a simple matter to show that the above commutation relations are satisfied if we 
define 

i ab = - l -[T a ,T b ] (43) 

Moreover, the matrices thus defined do obey the SO (5) algebra 

[£ ab , i cd \ = i(6 ac e bd + 5 bd t c - S ad £ bc - 5 bc t d ) (44) 
Let us adopt the following representation for the Clifford algebra: 

(45a) 
(45b) 
(45c) 
(45d) 
(45e) 

Then the charge and S z matrices take the form 

,51 _ 1 (o 2 \ „ 23 _ 1 f 1 



r 1 


= i <g> 


CT3 


r 2 


= (J i c 


3> o 2 


r 3 


= cr 2 C 


3> o 2 


r 4 


= cr.3 < 




r 5 


= -1 


® (J 1 



«"=2(o ;j <3 =Ho -ij (46) 

With the above generators Q and S z , the factorisation of the chiral field h in terms of the SU(2)2 P ' field g and of the 
charge boson <j) c must be 



10 



h = (,g T cos(a/7T^ c ), 9] sin(v90 c ), g t cos{^(j) c ), g± sm(y/Tv^) c ) (47) 
h = (,g T cos(Vtt0 c ),5t sin( v / 7T0 c ), 51 cos(V7r0 c ), ~5l sin(\/7r0 c ) (48) 

The explicit expression for the order parameter n a — tr(T a hh) in terms of the spin matrix field ( p8| ) and of the charge 
boson is then 

ni = Tr (G) cos(V^) (49a) 

n 2 = iTr (Gcri)sin(V^$) (49b) 

n 3 = iTr (Go- 2 )sin(V7r$) (49c) 

n A = iTi (GCT 3 )sin( % AF$) (49d) 

n 5 = - Tr (G)sin(0F6») (49e) 

We first notice that 721,5 form a spin singlet and are the real and imaginary parts of the complex d-wave order 
parameter Tr (G) exp(— iy/w9), whereas ^.2,3,4 form a spin triplet. We also see how the scaling dimensions of n^g 
diverge from those of n.2.3,4 when K c is different from unity: the fields cos(y / 7r#) and sin(-y/7T0) have conformal 
dimensions A = A = 1/(8K C ) while cos(y / 7r<I>) and sin(y / 7r$) have conformal dimensions A = A = K c /8. Thus 

A(n x ) = A(n 5 ) = 3/16 + l/8A c 

A(n 2 ) = A(n 3 ) = A(n 4 ) = 3/16 + A c /8 (50) 

We may also consider the SO(5) singlet hjii, which becomes simply Tr cos(y / ¥<i>) in this representation. This field 
is conjectured to be the charge density wave (CDW) order parameterBZI Within this model it has the same scaling 
dimension as the staggered magnetisation (or spin density wave), but in a more general model with a spin gap the 
two fields could have different correlation lengths since spin singlet and triplet states would not necessarily have the 
same excitation energy. 

To summarize, the components of the SO (5) order parameter have power- law correlations governed by the above 
conformal dimensions. When full SO(5) symmetry is present, A(n) = 5/16; and n is the vector primary field of the 
SO (5)i WZW model. When K c ^ 1, the SO (5) symmetry is broken and the staggered magnetisation is less (K c > 1) 
or more (K c < 1) relevant than the d-wave order parameter. The behaviour of the other possible fermion bilinears 
can be checked by a combination of Abelian bosonisation and an Ising model representation of bosonic exponents 
- we find that their correlations decay exponentially as a result of the gap in the flavour sector (for details see the 
appendix). Of course, since we are in one dimension, there are no real phase transitions, just enhanced fluctuations. 
Thus the fluctuations in the superspin channel are enhanced whereas other tendencies are suppressed. If a weak inter- 
ladder coupling were added to form a two- or three-dimensional system, then a mean-field treatment would lead to a 
phase transition in the channel that has the highest susceptibility (i.e., the most fluctuations), i.e. an ordered phase 
for the most relevant operator. Thus, this approach predicts d-wave superconductivity for weakly coupled ladders 
with attractive effective interactions (K c > 1). Finally, since the charge and spin sectors of this model are described 
by conformal fi^ld theories, one can also recover the finite-temperature behaviour of the correlation functions in the 
standard way.ES 

C. Lowest-energy states 

In the absence of interactions and interchain coupling, the low-energy sector of the model is especially simple: 
the theory is a SO (8)1 WZW model. The states fall into two representations of the Kac- Moody algebra: that of the 
identity which contains states with even charge, and that of the spinor (electron) field R^, which contain states of odd 
charge. Remember, this is just a complicated way of representing noninteracting fermion excited states. As SO(8) is 
broken into SO(5)i(8>SU(2)2', these representations break into a finite number of representations of SO(5)i(8>SU(2)f', 
as indicated in Eq. (|l^) . When the flavour sector is gapped, the low lying states must all be flavour singlets and so 
many of those representations become irrelevant, in particular all the representations of odd charge coming from the 
spinor R a . 

The only surviving Kac-Moody representation in the SO(5)i theory is that of the identity. Such a representation 
contains an infinite number of energy levels, and at each level the states fall into SO(5) multiplets. In the pure 
WZW model (before spin-charge separation) the excited states may be obtained from the vacuum by applying ladder 
operators associated with the SO(5) currents. Let us explain: in a system of finite length L, the currents may be 
Fourier expanded as follows: 
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l ab (x)=^2e 2ninx / L lf (51) 

n 

where the sum runs over all integers (positive and negative). From the commutation relations of the currents, one 
may infer commutation relations for the modes l® h and show that, for n < 0, Z° b is a raising operator for the energy in 
the WZW model. Of course, the Zg fc are nothing but the SO(5) generators and allow us to navigate within a multiplet. 

The multiplet content at each energy level may be easily obtained from the representation theory of Kac-Moody 
algebras, in particular by the method of affine characters.Ei Schematically, in the case at hand, the multiplet content 
may be expressed in terms of a spectrum-generating function X{q): 

X(q) = 1 + q 10 + q 2 {14 + 10 + 5 + 1) + g 3 (35 + 14 + 3 • 10 + 5 + 1) + • • • (52) 

where the coefficient of q A indicates the multiplet content of states with conformal dimension A. For instance, a 
term like 2q 3 ■ 10 in X(q) means that the multiplet 10 of SO (5) occurs twice with conformal dimension A = 3 in the 
right-moving sector. The full low-energy Hilbert space is a left-right product, encapsulated in the generating function 
X(q)X(q). For instance, the term Nq A 10®q A 5 stand for a left-right tensor product of multiplets, occurring N times 
at the energy level (2ttvf/ L)(A + A), with momentum (27r/L)(A — A) (vp is the common spin and charge velocity 
before spin-charge separation). 

The eigenvalues of S z and Q and the energy of each state in the right-moving sector may be encoded in a more 
general spectrum-generating function X(q, x, y): 



X(q,x,y) = Q A x 2S *y Q (53) 

states 



The advantage of spectrum-generating functions is that tensor products translate into ordinary products of functions, 
and direct sums into ordinary sums. Anticipating spin-charge separation, it is possible to write the function X{q.x.y) 
as a combination of spin-charge products: 

X(q,x,y) = X^(q,x)X^(q,y) + X^(q,x)X^(q,y) (54) 

where Xgp(q,x) is the spectrum-generating function for the spin-j Kac-Moody representation of SU(2)2 and Ac ' 2 ' 
is the analog for the charge sector. The lowest terms of these functions are 

xW>(q, x) = 1 + q(l + x 2 + x- 2 ) + q 2 (3 + 2x 2 + 2x~ 2 + x 4 + aT 4 ) +■■■ 

X$(q, x) = q 1 ' 2 ^ +x 2 + x- 2 ) + q 3 ' 2 {2 + x 2 + x~ 2 ) + q b ' 2 (A + 3x 2 + 3x~ 2 + x 4 + aT 4 ) + • • • 
Xi°\q, y) = l + q + q 2 (2 + y 4 + y~ 4 ) + q 3 (3 + y 4 + y' 4 ) + ■■■ 

Xi 2 \q, y) - q 1/2 (y 2 + y~ 2 ) + q 3/2 (y 2 + jT 2 ) + <? 5/ W + 2y- 2 ) + ■■■ (55) 

Again, the exponent of x is twice the value of S z and that of y is the charge Q. A term like Ax 2 y~ 2 q 3 in a j5^ ) 
would stand for four states with A = 3, S z — 1 and Q = — 2. The charge states represented in Ac ' have charge 
Q = (modulo 4) and those in Ac 2 ' have charge Q — 2 (modulo 4). From the above expressions and relation (S3), 
the full spectrum of energies and quantum numbers may be recovered. Of course, we must consider the left-right 
product X(q, x, y)X(q, x, y). That the expression (Q) is a sum of products, instead of being a simple product of spin 
and charge factors, means that one cannot consider the charge an d sp in spectra independently: there are "glueing 
conditions" between charge and spin states, conditions encoded in (p4j). 

When spin and charge separate, the energy levels shift in two ways. First, because of different spin and charge 
velocities, Xgp (q, x ) and Ac (<?, y) become respectively A^p (q s , x) and Ac \q c ,y), where q s — q v =/ VF andg c = q v "/ VF . 
Second, anomalous charge exponents change the conformal dimensions in the charged sector, whose structure deserves 
a more detailed explanation: excited states in the charge sector may be obtained either (i) by applying the creation 
operators associated with the charge boson <fi c (this does not change the charge Q) or (ii) by applying exponentials 
exp(z(5y / 7T0 c ) on the vacuum, where Q is the charge thus created. The generating functions in the charge sector may 
be expressed as 

XW(q,y)=X hos .{q) £ g^V 3 (56) 

Q=4r+l 

where r runs over all integers, I = or 2, and Abos.(<z) is the spectrum generating function associated with the boson 
creation operators only: 
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-^bos. (?) 



oo 1 



= 1 



2q 2 + 3q 3 + 5q i 



7q 5 



(57) 



When K c changes from its initial value of unity, left and right boson creation operators mix through some Bogoliubov 
transformation and the conformal dimensions associated with the exponentials of <f> c and <p c become 



A(Q,Q) 



1 


1 


32 




1 


1 


32 





(Q + Q) + VKc(Q - Q) 
(Q + Q)- Vk~c(Q - Q) 



Left-right products of spectrum-generating functions in the charge sector then become 

X^( qc ,y)XP(q c ,y)=X bos .( qc )X hos .(qc) £ q HQ,Q)qHQ,Q) y c 

Q=4r+i 
Q=4r'+£' 



(58) 



(59) 



The above expression, combined with Eqs (|54| , p5| , p8|) allows us to extract the energy, momentum, charge and spin of 
the whole low-energy sector for arbitrary values of u s , v c and K c . 

Let us consider, for instance, the first excited states. According to Eq. (52), They fall into the multiplet 10 of 
SO (5). The spin and charge content of such a multiplet is easily read from the corresponding S z — Q diagram of 
Fig[j]. The multiplet 10 consists of three spin triplets (of charge —2, and 2, respectively), plus a neutral spin singlet. 
After spin-charge separation and if K c ^ 1, the contribution of this multiplet to the spectrum-generating function is, 
according to Eqs (]5^,p5|), 



q c + 5,(1 + x 2 + x- 2 ) + 9yV (2 ' 0) g# (2 ' 0) (l + x 2 + x- 2 )(y 2 + y' 2 ) 
Thus, the energies of these states split in the following fashion: 

27T1V 



(60) 



E(Q = 0,j = 0) = 
E{Q = 0,j = l) = 
E(Q = ±2,j = l) = 



L 

2ttv s 

TTVs 
17 



~2L 



Kr 



1 



(61) 



The last of these states are in fact created by applying it operators (see|Eo. (34)). The energy levels are proportional 
to the scaling dimensions of the operators in the conformal field theoryBa 

To conclude, the eigenstates, in particular the lowest-energy states, fall into deformed SO (5) multiplets. The amount 
of deformation is exactly determined by the renormalised charge velocity v c and anomalous charge exponent K c . 



VI. THE LUTTINGER LIQUID CASE 



As we mentioned earlier, the case of two coupled Luttinger liquids is different to that of two Luttinger models. Let 
us consider, as an example of a Luttinger liquid, the one-chain Hubbard model at weak coupling U <C t: 

#Hub = ~t £ (4,a c r-+l,Q + 4+1,0^,0! J + U £ 'V.t'V.i (62) 
r,a r 

If we linearise about the right and left Fermi points as in (Q), and use the charge currents (||) and the corresponding 
spin currents 

a,P a/3 

we find the Hamiltonian density (vf ~ tao): 
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H H ub « -iv F J2( R t d * R ° - L\ t d x L a ) + j(f R +jt) + ^JrJl - - 2Uj R j L (64) 

a 

This Hamiltonian is not equivalent to a Tomonaga Luttinger model because the last two terms are not pure density- 
density interactions. The (j R + ] 2 L ) term will only renormalise the spin velocity v s and is not very important. 
However, the marginally irrelevant • jx, term cannot simply be absorbed in this way. It is this term which gives rise 
to logarithmic corrections to the correlation functions in Luttinger liquids - although otherwise it does not drastically 
change their properties, hence the utility of the Luttinger liquid concept. 

As this weak coupling bosonisation suggests, the coupling constant of this term is of the same order as the other 
couplings in the theory, so it cannot necessarily be jettisoned when we consider more complicated models such as 
the two-chain ladder. Let us therefore consider as a model for the generic Luttinger liquid two-chain ladder, the 
Hamiltonian (|l|) perturbed by marginally irrelevant spin current interactions in each chain: 

Wliq = H + TL mare 

^marg = ~A (jfil • hi +jii2 ' 3 Li) (65) 

where A > 0, and SriiSli are the right and left moving spin currents in chains £ = 1,2. It is instructive to write the 
perturbation in terms of the Majorana (vector) fermions £j. We find: 



j=l,2,3 



The first term is a marginally irrelevant interaction in the total spin sector (J/? = }m + jn2 as defined in (|6|)). But 
it is the second term which is most significant. It couples the fermions of the spin sector (£f = rightmoving, £f= 
lcftmoving, £ = 1,2,3) to one of the fermions in the flavour sector. So the spin and flavour sectors are no longer 
genuinely decoupled. 

Suppose that the flavour sector becomes gapped. Then: 

(UU) * (67) 

To first approximation, we can then replace i n ( |66| ) by its expectation value, and we see that the effect of a 
gap in the flavour sector is to generate a mass term for the fermions of the spin sector -j-a spin gap. This is a crude 
argument but it is borne out by the RG analysis of Refs. [l7| as well as numerical worked (in the notation of Ref. [l?] 

a finite backscattering corresponds to 7^ 0) which shows the existence of strong coupling regimes with a spin gap 
in a model of two Hubbard chains coupled by a small hopping. In general it is hard to estimate the size of this gap. 
If it is large, the low energy physics will be as described in Refs. [l7], Bll It could be, however, that in some models 
(with small A for example) the spin gap is very small, in which case for intermediate energy scales the behaviour will 
still be described approximately by the model dll). 



VII. CONCLUSIONS 



In this paper we have studied a system of two TL models coupled by a small interchain hopping. We have shown that 
this critical (at T=0) theory can be represented much more symmetrically than in the standard Abelian bosonisation 
representation as an SO(5)i WZW model, or equivalently as a system of 5 Majorana fermions, perturbed by symmetry 
breaking interactions. We have obtained the correlations of fermion bilinears in this theory and demonstrated that 
the components of the "superspin"EJ have power law correlations, and are enhanced, whilst other tendencies are 
suppressed. Conformal field theory allows us to obtain the exact energy levels in a finite size system and observe 
how the degeneracy of the SO(5) multiplets is broken by spin-charge separation (v c ^ v s ) and the presence of an 
anomalous exponent (K c ^ 1) in the charge sector. Except in the trivial noninteracting case, there is no exact SO(5) 
symmetry. In section VI we showed briefly how the inclusion of backscattering results in the appearance of a spin 

g a P- n 

In the light of these results, S.C. Zhang has recently shownEJ that a whole class of ladder systems with more general 
interactions have Hamiltonians with microscopic SO (5) symmetry. Knowing the Hamiltonian does not tell us about 
the strong coupling behaviour at low energies. But a continuous symmetry like SO(5) cannot be spontaneously broken 
in one dimension and must therefore be present in the low-energy theory as well. The latter must be described by a 
SO (5) WZW model, perturbed by various primary fields, perhaps with a critical point or line in the space of coupling 
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TABLE I. Fermion Bilinears. 



constants (a conformal field theory with Lie-group symmetry is necessarily a WZW model) . This is the simplest class 
of low energy theories with SO(5) symmetry in 1+1 dimensions. 

The SO(5) symmetric description of ladder models, which are clearly related to popular, models of the Cuprates, 
and the similarity of the form of the theory in ladders to that proposed by S.C. ZhangEl for the 2d Cuprates, is 
certainly encouraging and suggestive. Nonetheless, in view of the many special features of one-dimensional theories 
we are cautious about drawing more general conclusions. 
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IX. APPENDIX: ABELIAN BOSONISATION OF FERMION BILINEARS 

We want to know the long distance (low energy) behaviour of the correlations of various bilinears. In Sect. (|y|) 
we demonstrated that the correlations of the "superspin" order parameter could be deduced within the framework of 
non-Abelian bosonisation from a careful analysis of the operators of the conformal field theory. We can further justify 
this analysis and find the behaviour of the other fermion bilinears explicitly by using Abelian bosonisation. 

We introduce an Abelian boson <$> l a for each species of fermion i = 1, 2, a =|, 4- Then we introduce linear charge 
and spin, bonding and antibonding combinations of these fields: 

K,s = ^m±n) (68) 

If one carefully applies Abelian bosonisation to the original Hamiltonian ([!]) , taking full account of the anticommutation 
factors£il one can identify each of these Bose fields with two of the Majorana fermions introduced in section (4). The 
field is simply associated with the two charge fermions, the $J field represents two of the flavour fermions. $+ 
represents two of the spin fermions, and $J comprises one flavour fermion and one spin fermion. 
After bosonisation, the various fermion bilinears have the form 

O - e - l ^ x (69) 

where the different Xs are given in Table J. We will briefly describe how we arrive at the long distance behaviour of 
their correlations. We find straightforwardly: 



A(e ±!V ^) = A(e ±lV¥e = + ) = - 



A(e ± ^*^) = — 
8 

A(e ±I ^ + ) = _L (70) 

(The scaling dimension is D = A + A and here A = A so D — 2A) But the charge and spin (— ) fields are a little 
more subtle. In our model, the flavour sector acquires a gap. Since we start with K c < 1 this strong coupling regime 



15 



corresponds to the limit K c — > in the c— sector, whence to leading approximation we can replace the complex 
exponents by their expectation values: 

(e i ^") = (71) 

Higher order corrections will die away exponentially. Thus the bilinears (a) , (e) and (f ) in Table [| die away exponen- 
tially and the corresponding tendencies are suppressed. 

The exponents of $J are a little more subtle since we know that only one of the Majorana fermions to which <6~ 
corresponds is gapped, whilst the other remains gapless. Here, , howe ver, we can make use of their representation in 
terms of the corresponding Ising order and disorder operators-r^rTp 5 ! Introducing Ising order and disorder operators 
<Jf,iif corresponding to the Majorana flavour fermion £3 and er s ,/i s corresponding to the spin fermion we can 
identify the following approximate operator correspondences: 

COSy/n^J ~ a f a s 

sin \Ztv$J ~ fj.fH s 

COS ~ CTfUs 

sin -^#7 ~ M/^s (72) 

When the flavour fermion becomes gapped, either (it/) = 0, (jii/) ^ or (erf) ^ 0, (///) = 0, depending upon the 
definitions. Thus to first approximation these can again be replaced by their expectation values. The Ising model 
corresponding to the spin fermion £^ remains critical and has scaling dimensions A = A(cr s , fx s ) — 1/16. In this way 
we obtain the long distance asymptotics shown in Table || The only bilinears which still have power law correlations 
are the interchain pairing fl39|), represented by (b) in Table |, and the staggered magnetisation (|37|), represented by 
(c) and (d). Thus, it is precisely the components of the unified order parameter n a that have power law correlations, 
while all the other tendencies around ±2fcp are suppressed. Note that the scaling dimensions that we find agree with 
those found in Sect. from non-Abelian bosonisation (cf. Eq. (p0|)). 
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